Square patterns in Rayleigh-Benard convection with rotation about a vertical axis 
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We present experimental results for Rayleigh-Benard convection with rotation about a vertical 
axis at dimensionless rotation rates < Q < 250 and e = AT/AT C — 1 < 0.2. Critical Rayleigh 
numbers and wavenumbers agree with predictions of linear stability analysis. For £1 > 70 and small 
e the patterns are cellular with local four-fold coordination and differ from the theoretically expected 
Kiippers-Lortz-unstable state. Stable as well as intermittent defect-free square lattices exist over 
certain parameter ranges. Over other ranges defects dynamically disrupt the lattice but cellular flow 
and local four- fold coordination is maintained. 
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PACS numbers: 47.54.+r,47.20.Lz,47.27.Te 

The elucidation of spatio-temporal chaos (STC) re- 
mains one of the major tasks in the study of pattern 
formation in nonlinear dissipative systems. [|l| The best 
opportunities for theoretical understanding of experi- 
mental observations of a chaotic state exist when the 
mean-square amplitude of STC evolves continuously (i.e. 
via a supercritical bifurcation) from a spatially-uniform 
ground state; for in that case it should be possible at 
least in principle to derive from the equations of motion 
of the system a systematic, simplified description in the 
form of Ginzburg-Landau equations. However, exper- 
imentally accessible supercritical bifurcations from the 
uniform state to STC are rare because most systems with 
supercritical primary bifurcations become variational as 
the threshold is approached from above and thus ap- 
proach a time-independent state. Convection in a shallow 
horizontal layer of a fluid heated from below (Rayleigh- 
Benard convection or RBC) and rotated about a vertical 
axis is one of the exceptions. Convection occurs when 
the temperature difference AT exceeds a critical value 
AT c (0) (Cl is the rotation frequency) and leads to a ve- 
locity field v. The Coriolis force Oxv renders the system 
non- variational even close to the convective threshold and 
thus permits the existence of STC at onset. Kiippers and 
Lortz ||] (KL) predicted that a primary supercritical bi- 
furcation leads to a state of unstable convection rolls pro- 
vided Q > Qkl- In this KL state, rolls of one orientation 
are unstable with respect to another set of rolls with an 
angular orientation relative to the first which is advanced 
in the direction of rotation by an angle 6. ||,||,f| The new 
set, however, is equally unstable to yet another, and so 
forth. Several experiments confirmed the existence of the 
KL state at relatively krge e = AT/AT C - 1. Re- 
have shown its existence close 
12 < fl < 20 and have verified 



cent experiments 10 
to threshold for Qkl — 
the supercritical nature of the bifurcation. The instabil- 
ity was observed to lead to a chaotically time dependent 
co-existence of domains with different roll orientations. 
Two examples [|l0[ of this domain chaos are shown in 
Figs. [l|a and b. 




FIG. 1. Shadowgraph images of convection patterns viewed 
from above for CO2 (a and b, only parts of patterns are shown; 
adapted from Ref. Jl(]]), water (c and d, the entire cell is 
shown), and Argon (e and f, 90% of the cell radius is shown) 
with (a, H, T) equal to (0.93, 19.8, 40), (5.4, 170, 4.8), and 
(0.69, 181, 8.3) respectively. For small Q (a,b) the patterns 
have domains which are typical of the KL instability for all 
e near zero, as illustrated for e = 0.06 (a) and 0.18 (b). At 
larger Q, square patterns occur close to onset [(c) e = 0.09 
and (e) e = 0.04], but states similar to the KL domains are 
observed for e > 0.1 [e = 0.12 in (d) and e = 0.13 in (f)]. In c 
and d, the cells along the periphery are the wall mode. 



Theoretically it is expected P,p|,^,pT[ that the nature of 
the bifurcation and of the nonlinear state above it should 
remain qualitatively unchanged as Q is increased above 
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the range explored previously by experiment. The mea- 
surements which we report here show that this is not the 
case. For ft = 2nfd 2 /v > 70 (/ is the rotation frequency 
in Hz and v the kinematic viscosity) we find that the bi- 
furcation does remain supercritical as predicted, but that 
the convection pattern above onset has no similarity to 
the expected KL state. Instead the pattern consists of 
cells which are usually arranged so as to have local four- 
fold coordination. Over significant parameter ranges the 
cells "crystallize" and form a square lattice. Q Typi- 
cal examples are shown in Figs. |l|c and e. Depending on 
parameter values, the lattice can be stable, can be inter- 
mittently disrupted by defects which seem to be injected 
from the boundary, or can be continuously disorderd with 
many defects within it (maintaining, however, the cellu- 
lar character with predominantly local four-fold coordi- 
nation). At larger e (e > 0.1), we do observe patterns 
reminiscent of the KL state, as illustrated in Figs. 0d 
and f. The existence of the square state for e-values be- 
low those where the KL state is observed is contrary to 
theoretical prediction. This disagreement with theory 
occurs in a parameter range where the stability analysis 
was thought to be reliable, and thus presents a significant 
unresolved problem in the field of pattern formation. It 
also implies that the use of the KL state as a prototype 
for STC at onset is valid only at relatively small Vt. 

We used water and Argon in two different apparatus 
]l3| , p!4| . The Prandtl numbers a = vj k [k is the thermal 
diffusivity) were 5.4 and 0.69 respectively. The tempera- 
tures of the tops of the cells were regulated to ±lmK at 
30.07 ° C (water) and at 37.5 ° C (Argon). Aspect ratios 
r = r/d = 4.8 and 8.3 (r is the radius and d the height 
of the fluid layer) were used for water (r = 38.1 mm, 
d = 7.9 mm) and Argon (r = 33.25 mm, d = 4.0 mm) 
respectively. The apparatus were rigidly mounted on ro- 
tating tables capable of rotation up to / = 1 Hz, cov- 
ering the range < < 250. Shadowgraph assemblies 
were mounted ontop of the convection apparatus to ob- 
tain images in the rotating frame. For Argon, the pres- 
sure was controlled to ±2mbar. The experiments were 
carried out by raising the bottom-plate temperatures at 
fixed Q. Both heat-transport measurements and shad- 
owgraph images were taken after waiting for at least a 
horizontal diffusion time Th = T 2 d 2 j k. Onset was deter- 
mined both from the heat transport and from amplitudes 
of the shadowgraph images. 

The dimensionless control parameters are the Rayleigh 
number R = agd? /ST / nv and Q. Here a is the isobaric 
thermal expansion coefficient, g the acceleration due to 
gravity, and AT the temperature difference across the 
cell. The effect of centrifugal force given by the ratio 
of centrifugal-to-gravitational force F = 47r 2 f 2 r/g is less 
than 0.03 in most of the experiments and maximally 0.12 
for a few runs. Since Q is scaled by r v = d 2 jv, increasing 
the pressure P of Argon increased ft because it decreased 
v. We used this to vary Q by almost a factor of 4 at 



fixed / (and thus F). Increasing P changed a only from 
0.68 at 20 bar to 0.695 at 40 bar. It decreased the pa- 
rameter Q jL5| which describes the extent of departures 
from the Boussinesq approximation. For Argon Q ranged 
from 0.03 at 20 bar to 0.008 at 40 bar for the highest Q 
(largest AT). For water Q ~ —0.02 at large Q. Since Q 
was always small, we expect non-Boussinesq effects to be 
negligible. 
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FIG. 2. The critical Rayleigh numbers R c (Q)/R c (0) (top) 
and wavenumbers k c (Q) / k c (0) (bottom) as a function of Q for 
water (a = 5.4, open circles) and Argon at 20 bar (a = 0.68, 
open triangles), 30 bar (cr = 0.69 inverted solid triangles), and 
40 bar (cr = 0.695, solid triangles). The results agree with the 
predictions for a laterally infinite system (solid lines) Jl(| . 

Figure || summarizes measurements of the critical 
Rayleigh numbers R C (Q) and wavenumbers k c (Q) for the 
onset of convection in the sample interiors (the "bulk 
mode"). The results agree with predictions based on lin- 
ear stability analysis for a laterally infinite system. Jl(| 
These predictions are independent of a and are given by 
the solid lines in Figs. ||. For Q > 70, the bulk mode is 
preceded by a "wall mode" consisting of a wave traveling 
(in the rotating frame) in the direction opposite to Q. 

The wall mode persists above R c , and can be seen in 
Figs. 0c and d. The primary object of the present paper 
is the bulk mode. 

Figure || shows shadowgraph images for e = 0.04 for 
various Q in Argon. At small Q (Fig. ^|a) the contrast of 
the image is poor because AT and the wavenumber were 
small. Nonetheless one can see that almost the entire 
cell was occupied by a single domain of rolls. However, 
the time dependence of the pattern clearly revealed a dy- 
namics characteristic of the KL state. PJnl At somewhat 
larger Q (Fig. ||b) several KL domains exist simultane- 
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ously in the cell. These observations are consistent with 
previous work |9],|l0| using larger T and Q < 20. This 
situation changes when Q is increased beyond about 70, 
as illustrated by Figs. ||c and d. For Q = 73, the pat- 
tern is made up of cells rather than rolls. There are 
large domains where the pattern has square symmetry. 
Six domains of squares meet in the cell center and locally 
enforce a six-fold symmetry, maintaining however the cel- 
lular character of the flow. At larger Q — 145 the pattern 
consists of a square lattice, with some domain walls and 
defects (right part of Fig. ||d) apparently provoked by the 
cell wall. 




FIG. 3. Patterns for Argon and various O at e ~ 0.04. Q 
values are (a) 13, (b) 40 (c) 73 and (d) 145. Only 90% of the 
cell radius is shown, c and d are at the same physical rotation 
rate and at different pressures. 




FIG. 4. Disordered cellular patterns for a = 5.4, Q, = 170, 
and e — 0.052. The entire cell is shown. Along its periphery 
is the wall mode. 



There are ranges of e, and a over which the patterns 
near onset were dynamically disrupted by many defects, 
although cellular flow, often with local four-fold coor- 
dination, was maintained. Two examples for a = 5.4, 
Q = 170, and e = 0.052 are shown in Fig. ||. In order 
to study the pattern dynamics in greater detail, we ex- 
amined the time evolution of the structure function 5(k) 



(the square of the modulus of the Fourier transform of 
the pattern). Two examples of *S(k) are shown in Fig. ||. 




FIG. 5. Structure functions of images lc (right) and Id 
(left). 
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FIG. 6. Angle- time plot for the radial average of the struc- 
ture function for a = 5.4 at Q, = 170 for (a) e = 0.072 and (b) 
e = 0.091. (a) is a representative section of a long run. The 
system intermittently changes between a perfect square lattice 
and a disordered pattern with local four- fold coordination. 
For (b) e was stepped from 0.14 where the pattern was disor- 
dered (KL like) to 0.091 at time t = — 3 hours. After about 
16u, the system spontaneously chose the perfect-square pat- 
tern which then persisted for more than 38r^. 

The left one is typical of the KL state, whereas the right 
one represents a square lattice. For Q = 170 and a = 5.4, 
the radial average 5(6) of 5(k) is shown in Fig. | as a 
function of time. Figure ||a is a section of a long run 
at e = 0.072; we believe it is representative of a static- 
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tically stationary process. Disordered regions alternate 
irregularly with well ordered square lattices. We believe 
the disorder is provoked by defects which are injected by 
interaction with the wall mode, but further work is re- 
quired to examine the mechanism. Figure ^|b is the result 
of an experiment where the system was kept in the KL 
regime at e = 0.14 for a long time, and where at time 
t = — 3 hours e was reduced to 0.091. The system re- 
mained KL-like for nearly 44 hours (about 16 horizontal 
thermal diffusion times), and then suddenly crystalized 
into a square lattice which remained stable thereafter. 



impurities or poorly conducting boundaries cannot be in- 
voked to explain our observations. Thus the occurrence 
of square patterns in rotating RBC at small e is unex- 
plained, and we conclude that a real discrepancy exists 
between theory anc ^ experiment. 

We are grateful for stimulating discussions with F.H. 
Busse and W. Pesch. This work was supported by the 
U.S. Department of Energy Grant DE-FG03-87ER13738. 
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FIG. 7. The rotation rate u) of the square lattice for 
Q = 170 and a — 5.4. Time is scaled by d 2 /v. 

It is interesting to note that the square lattice, when- 
ever it exists, rotates at a rate co relative to the con- 
vection cell in the direction of the overall rotation SI, 
albeit extremely slowly. Results for uj are given in Fig. 0. 
The measurements suggest that uo vanishes as e goes to 
zero, as indicated by the dashed line. There is of course 
no theoretical prediction for a;, since even the existence 
of the lattice is not understood. For comparison, under 
similar conditions the wall mode is travelling in the op- 
posite direction at the much faster angular velocity of 
about — 0.15z^/<i 2 , nearly independent of e. We believe 
that the interaction between the fast-moving wall mode 
and the slow-moving bulk mode leads to the defect for- 
mation which tends to disrupt the square lattice for some 
parameter ranges. 

Stable square patterns occur under other conditions in 
RBC. Convection of a pure fluid in a cell with insulating 
top and bottom plates jl^] and binary-fluid convection 
with small Lewis numbers and a positive separation ra- 
tio ]l9| are two such systems. However, in those cases the 
critical wavenumber is reduced. Our critical wavenum- 
bers agree with the linear stability analysis for rotating 
RBC of a pure fluid with perfectly conducting top and 
bottom boundaries. Our top and bottom boundaries 
have conductivities which are orders of magnitude larger 
than those of the fluids. Our Argon samples, even if 
they were contaminated with another component, would 
have Lewis numbers of order one |2(]] and thus would not 
produce a square pattern due to mixture effects. Clearly 
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